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Many- Particle Theory of Nuclear Systems with 
Application to Neutron Star Matter 

This is a report on the research carried out on the followijig 
projects since submitting the semi-annual status report on 
April 5, 1973. 

A. Calculation of an improved energy- density relation 
for the normal state of neutron-star matter 

B. Calculation of the effects of superfluidity and 
polarization in neutron star matter 



project a 


The Calculation of an Improved Energy-Density Relation in Ncutroa 
Matter in the High- Density Region 

A-1 Constraints on Variation 

The theoretical formalism is outlined in section III.l of the orig!'-'' ' 
proposal. In particular, the energy of the normal state of neutron matce.>, . 

£ , is developed in a cluster series 

IV 

= Sp + ^2^ £3 • 

In this report, we shall use €,^to refer to our approximation for energy, 
namely, ( i- £3 ) . 

Reference is made in the original proposal to certain physically 
motivated necessary conditions on the radial distribution function. These 
in turn, constrain the variational parameters in the trial two-body 
correlation factor f(r). There is also the so-called Pauli condition which 
restricts f(r) directly. It arises from the effect of the Pauli Principle 
which prevents particles in the fermi sea from scattering back to the 
occupied states , The conditions which we use as constraints on our 
variational procedure are listed below: 


(I) 

H 

H 

0 

(Pauli Condition) 

(II) 

SCVO^o 


(III) 

Xt ^ »-243 

(Coulomb Inequality) 

(IV) 

so 


(V) 

5Co> = 0 

(Srructure-Factor Simi Rule) 


I^ and £^^are defined as follows; 


where Ig, S(K), 
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£(x) - ^ r S^/»<.x CoS X N 
X —2, — ) 

SCk.) - l + •+ S^‘Vi«.) H' • •■ 

S^"tk:> = P J <ij: \ ) 

S^\k-> ^ 

^Iv-) =: f\r) ( ' “ \ 

0(9 IP^S 

S^"V) - ~(i^)TTf f£r> Jas [dt St r^cs:. . 

^ 0 ir-s( 

• iCkpW Icfefii ICkpt)^ 

1. «• r+ s 

+ TTf iS^ Cjls 'Vjvs> T^tt) • 

- c iW» + 

(OT T* ^ 

Xc ■» Ic + ^c -*■'•• 

= (;^^Ttp■)^^ (dl'T T ( ' - ' 2 [“'>'>) 

A- r s{cr-^ 

^ c) 

iijr = ■£ f |4r ( * "‘i t^ckpr)) 
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We note that 

Z (AfiO = i + s ‘“(w 

Thus condition CiV) is equivalent to S^°^(0) = -1. Condition (V) means 

f 

S(«> = 1 + S®(0) t + = 0 . 

we use it in the truncated form 

J + -V =0 

The origins of conditions (II) and (III) are discussed in E, Feenberg, 
Theory of Quantum Fluids (Academic Press, New York, 1969). 

The relative order of magnitude of the terms in the cluster series 
for energy and all other associated cluster expansions is determined by the 
"correlation parameter" ^ defined by 

J - nur , uj ^ J^Cf*(Y) "I ) 

We note also that the five conditions are not completely independent of 
each other because of the following relations: Condition (V) is related 

to the k = 0 version of (II), The former can be satisfied by first 
satisfying (IV) and then requiring also that S^^^(o) = 0, 

In order to expedite the manerical calculation, we usually avoid im- 
posing (II) and (III) directly. Often we find that (V) enables us to find 
more easily the region of the parameter space where (II) and (III) are 
satisfied. Similarly, (IV) (or (I)) enables us to locate regions where 151 - 


small , 
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A-2 Numerical Procedure and Results 

Several methods of calculation', each involving a different set ot 
criteria governing the choice of constraints in the variational procedure, 
have been pursued. As indicated in the proposal, we have used one- 
parameter, two-parameter and three-parameter forms for the trial two-body 
correlation factor f (r) . Calculations using the various methods of 
approximation have been carried out at the density corresponding to fermi 
wave number, = 3.5 fm'^ in order to test the reliability of these method; 

r 

in the high-density region. The two-nucleon potential used for the 
purpose is the hard-core potential of core radius 0.4 £m given by Ohmura, 
Morita and Yamada (omy-4) (Progr. Theoret. Phys. 1^ (1956) 222). This is 
a central potential of the form 


a , A = r. 


^3 “^4 ~ 


- 3 f Ch.o; 


V, (r:,) = 


oo 




v.i. e.'^* 


where i = 1, 2, 3, 4 denotes, respectively, the component appropriate to the 
singlet-odd, singlet-even, triplet-even, and triplet-odd state of the two- 
nucleon system. The A^are the corresponding projection operators. The 


parameters are 
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V 


o2 


= -235,41 MeV, ^ 


.0344 fm 


-1 


-1 


-13 


V , = -475.04 MeV, 3- = 2.5214 £m 
o3 -3 

''ol ‘ ''o4 ■“ ? 

y* =0.4 £m (in all States) 
c 

These parameters are chosen to £rt the £ollowing data characterizing tn.- 

low-energy interaction o£ two nucleons in free space: 

Binding energy of the deuteron = 2.226 MeV 

Triplet scattering length of the neutron-proton system = 5.378 x 10 

Singlet scattering length of the neutron -proton system 23.69 x 10 cm 

-13 

Singlet effective range of the neutron-proton system = 2.7 xlO cm 
A two-nucleon state must be either singlet-even or triplet-odd 
according to the Pauli Principle. Therefore, in our calculation, there is 

no need for 
Method 1 


CUI 


A trial two-body correlation factor. 


^Or) ■= J 


-.yj tr-rc> 


r 6 Tc. 

I - c ' , 

has been used. According to the variational principle, the minimimi of 
obtained using any trial wave function provides an upper bound to the true 
energy E. Therefore, at each density we minimize respect to 

yL/ . If this minimum occurs at ^ , then C /^) is calculated. The 
corresponding approximation to the eneirgy .per neutron is 
The results obtained by this method are given in Table A-1. 
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Method 2 

The following two-parameter form for £(r) has been used in this method 








IS e 




1 


Y- < 


j 


r > r*. 


The additional parameter gives £(r) more flexibility in order to help satisfy 
the constraints better. Several alternate approaches in imposing these 
conditions have been used. These are described below; 

(A) The parameter is fixed at 2.0 £m Y is varied and (fjp+^ 2 ,) is 

found to have a minimum with respect to IS at V = is 

calculated. 

(B) X is determined by condition (I). Then the minimization procedure of 
method 1 is attempted. 

(C) The same as CB) , using condition CTV) instead of condition (I) . 

(D) The same as (B) , but using condition (V) instead of (I), followed by an 
atten^t to use the minimization procedure of method 1. 

In methods (B) and (C) neither (Ep + £ 2 ^ ^2 ^ 

to have a minimum with respect to ^ . Results for method 2A, 2B and 2C 
presented in Table A-1 are intended to illustrate the following fact: As 
the value of fJ increases, energy gets larger and convergence of all cluster 
expansions improves. This causes some ambiguity regarding the determination 
of the minimum energy. However, the ambiguity can be removed when we impose 
the conditions (11) and (III) on the €^(//) curve. This is done in method 
2D for which the complete results are given in Table A-2 for ^ values in the 
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range 2.0£ 7.0 It should be pointed out here that there are 

two values of V that satisfies S(0) = 0 in method 2D. Only the smaller of 
these values of % is used for each , since that corresponds to lower 
energy and better convergence of all the cluster expansions. 

Method 3 


Here we adopt the following three-parameter form for f (r) ; 


•f(lr) 





X C 


- >i.i>p-rcV 




)r < n 
r >r«. 


Several different procedures involving this f(r;) have been attempted. Of 
these, the ones that turned out to be most fruitful are described below. 

(A) For given // and X , is determined by conditinn Cl) . Then it is 

found that at eachX, both and €^have a minimum with respect to 

y at approximately the same value Now we bring the results for 

^ Cp >10 better agreement with the constraints by changing Y to the lower 
of the two values that will make S(0) = o. is to be determined by a 
further minimization of energy. The results are presented in Table A-3. 

(B) In this method, the procedure in (A) is followed upto the point of 

obtaining the I (/^JO results. Then X is chosen as the value ^ which gives 
the smallest Hi. Then is recalculated and its minimum with 

respect to ^ is sought. The results are given in Table A-4. 

(C) In this procedure, again we let and X vary, but use condition (I) to 
fix ^ = y as a function of ^ and X . Then at each trial value of ^ , X is 

chosen as the valued? for which UXt is smallest. The resulting data for 
£ and other auxiliary quantities are shown in Table A-5. The lowest 
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energy subject to the constraints is to be determined from this data. 

(D) This procedure is the same as (C) except for the following changes; 
Condition (IV) is used instead of (I) to fix Y ; is chosen correspond- 
ing to the smallest value of 1 Ig\ . The results are given in Table A-6. 

A- 3 Discussion 

The lowest energy obtained by each method is shown in Table A- 7. 

The criteria that are used to determine the lowest energy are the following: 

(I) Whenever ^ (^) or l^as a minimum with respect to ^ , we 

choose that as the lowest value in spite of the violation of S(0) = o that 

is usually associated with it (but, only in the high density region). 

However, S(k) > 0 is satisfied for k values not near k = 0. This can be 

considered adequate because in the high-density system the larger k 

values are much more significant. Besides, in a convergent cluster 
expansion for S(k), 

S(k) = 1 + S^°^(k) + S^^^(k) + ... 

the neglected higher order terms, though expected to be small, may never- 
thless be sufficient to "repair'* small violations of S(o) =0 and 
S(k) i 0 for small k. In short, the results for energy obtained from a 
calculation in which all the cluster expansions converge rapidly and all 
the conditions are satisfied with the exception of S(k)>o at small k, 
may be considered reliable, (2) Even when (£p + 62) or no minimum 

with respect to /J in an unconstrained variation, it is found that conditions 

(II) and (III) cannot be satisfied when is below a particular value ^ , 
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Thus £ (|^) may be taken as the lowest energy consistent with the consrra;. r .s 
(3) Note that condition (II) is imposed in the modified sense discussed 
above, except for methods 2D and 3A, where we have sought to satisfy ^ 

for all k. The results form these latter two methods, when compared with 
results from other procedures, give us an estimate of the error we ma/ 
allowing through the violation of S(k)>o for small k. (4) The reliabil^i,^ 
of all our results depends on the size of obtained; if\?lis large, 
adequate convergence of the cluster development is in doubt. Based on 
all these criteria, so far the best results are obtained from method 3D. 

we are in the process of testing four additional procedures involving 
the three-parameter correlation factor. This is expected to^be completed 
in about three weeks from now. Then we will adopt the most reliable of 
the methods we have tested and carry out the complete calculation for 
the entire density range 0.2S£ kp£ 3.5 It should be enphasited 

here that the difficulties associated with imposing S(k) > 0 for all k 
do not arise in the intermediate - and low-density regions. 



MeTHOP 


MeTHOl 

, 5 

H 

cW) 

> 

S(o> 

A>r 


1 

■ -3-t>o 

^•3 

O 

+?-t 

-1 57 

-0-1? 

2A 

-l'3o 

I'O 

J 20 

-U| 

-0'4l 

-O'Ob 

2B 

0' Ilf 

10 

1-30 

-l-S 

0-30 

o 

26 

0- 10 

A'O 

‘•53 


o*i? 

0 

2C 

-O'H? 

2-0 

I'Zfc 

-20 

o 

-O'Oi 

2c 



l'40 

-0'^ 

o 

-O-oi 


^2. ^3 

(Mex> CM-fcV) CMeV) (fAtV) 

-bO»l 9i‘Z 6^-1 (bO. 

0'?t -S&'l 13-*5 [fol. 

1-31 “ 51-2 

\'Ol +\%‘l l^l. 

I'ZO “5t.2 

»' 01 > + 12-7 12>0 237 


2B and 2C 




t 

t^vwT') 

I 5 toi 

S Co> 

SCo) &K 



(Mev) 

2Tfr 

c 1 V\W) 

^3 

C|H.fcV) 

C»v\c\» j 

SCk'^^o . te>o 

-l%0 

20 

I.IU5 

1-37 

0 » 0 & — * 0 ‘ (p5 

-O’Oi 

OSfc 

-S9-3 

43.0 

72'$ 

t r 


-l«o? 

3-0 

1 ’ 2 of ^ yL 

0-b7 

0'O» -0-33 

-0 0^1 


-41*1 

III • 

75 '4 

l%7. 

5*0 5. k ~ T 

-0 74 

j 

4'0 

-'-^3 

0-^3 

0.00 -p.2| 

-|'D3 

O' 9 b 

-4'4I 

14$. 

*74'0 

224. 

M o vne 

-0-SI 


l'3U -1*31 

0 l\ 

000 -O-l? 

-D-*3 

® '94 

+9?.8 

l«lb. 

77 'S 

274. 

None 

-O'HT 

4 0 

\ -123 

0’Z3 

0*00 -0‘»Z 

-'P'OZ 

C>'9i 

103. 

as?. 

$10.1 

335. 

None. 


7^0 

-i-n 

0*17 

0*00 'O'Oi 


0'«9 

17?. 

111. 

^3.4 

4 n. 

l\l c?v> ^ 


TABLE A-^2 


Results from Method 2D 



1^ 

(f«') 

^ S^?o) 

SCO) Ajsf Ig 

Tc 

£z 

(««v) 

fp+Sj, 

tMev) 

crAtv) 

fv, 

CtA«V) 

S(it) BO , K>o 


I'O 

20 

OJM -31.-5<1 

i*5^ -0>O» -0*17 -O'O^ 

0*2? 


^2*7 

10 'k 

IfcJ* 

zijwr'z Wi ? 

-1*21 

^•0 

3-0 

0-52 -\’%0 

0'?l +o»o» -o*Mo - 0*05 

0'8«» 

-4q*5 

103* 

76? 


/j*0f^*f k ? 


5«0 

^'0 

0'6i 

0*<47 -O'OZ “0*24 “0‘0i| 

©••It 

-<g*2 

13^. 

TT'6 

212. 

Wo V^C 

-D‘4,3 

G ' 0 


0*«2 -'1-3+ 

0 34 O'OO -0'17 '0-03 


+ 26.1 

ITS. 


ZS7. 

Nov\c 


7'0 

6*0 

0'% ■*'l'29 

0'29 0*00 “0*13 -O'Di 

0-?2 

73*1 

;i35. 

^1*2 

317. 

Novit 

-O'ifl 

?'0 

70 

I-IO -I'H 

0*14 O'OO “O'lO -O'OZ 

0-sn 

15Z. 

3o4. 

84*7 

3?^. 

Nov»^ 


TABLE A- 3 


Results from Method 3A 



S"co^ S Co) 5 to; Ajr I’b ^2- Xk.):5^o 

tf^fV) CiAtV; Crv>ev) (fAeV) 5«»A»sfi:ei fo 


l\\ 

193. 

414 

Z^S. 

^ k < 

-H'3 


100 

2 \\. 

k < 

-2^.1 

\Z9. 

m-2 

202. 

1 t k < 

-21 0 

\29. 

19.1 

200 . 

t •& '< k. ^ 

“25't 

Ul . 

7S‘i 

20Z. 

ckc 

-2h«} 

130 - 

l^.S 

209 . 

i * G fvo^ ^ k < 


TABLE A- 4 

Results from Method 3B 





f ^ if Sfe> AK 

If"”*) 


^2. 

C<^e») CMev) 


•3 


SCk")^© 


Crnev) CMeV) 


0*23 

\’H 

1-75 

2 3t 

-0 3b -2-52 

“1*14 

0*32 

0 

Q.%\ ^i+5.^ 

lob' 

85*5 

0 *4 < k < l*2Siii 

142 • ^ i_ , -1 , 

2-0 < k<- 3*1 bx 

o.<s 

t' (« 

2-00 

2-33 

-0 J4<^ -2'00 

”I*HH 

0*28 

0 

.* 

l*o<1 -31-0 

l»3. 

6u 

80Z 

O'b < k< 1*2 5*w>* 

Iff*! 

2*0 4 k. ^ 5*0 

0'«7 

\>% 

22t 

242 

-0.M7 -l'b3 

-M2 

0-2$“ 

0 

UH ”2 9^0 

124, 

154 

0*7 <. k 4 1*3 b**' 
22<^<T ff' 

O-R 

Z-0 

2'€»0 

2-63 

“«43 H-37 

-040 

0-23 

o 

t'OH -127 

151* 

13*5 

ZZ$^ K>G'7^..\ 

07^ 

2'Z 

2-SO 

2'55 

H>-Sb 

“045 

0*22 

0 

1*12 348 

»5b. 

bl .4 

225. k>0*4fm‘ 

0 -17 

2*5 3'2Z 

2 *6(> 

-0'6o -0-9^ 

-041 

OZO 

o 

»'07 34*9 

187 


248. k>\'0^>-r‘ 

0‘(t 

Z'« 

3 '5? 

2-51 

”0*62 —0 ’S(o 

-6 49 

<5I<J 

0 

1*05 52*1 

2o4. 

564 

261. k > 1* 1 

0 -17 

3-0 

3 '36 

2'fel 

”045 

-0-HZ. 

0 13 

o 

<•05 7o*(i 

223. 

54.0 

217. k<>i.tfwr’ 

0-17 

3-2 


242 


-O’it 

0*19 

0 

I'OI 98’2 

2^1. 

s-ii 

242 . k > 1 . c» f ’ 

O'lS 

3'? 


273 

— 0 45 “itf'bS 

-0-33 

o*n 

o 

047 »28. 

281. 

4fc*3 

321. k>0'*‘?fW^ 

d-l$ 

3*% ^ 

272 

-0*6b -643 

-6 21 

0*17 

o 

6 15 /55- 

3o8, 

444 

352. k704fUT’ 







TABLE 

A- 5 






Results from Method 3C 




t 

y y 

^%'i SCo) 

5ro; 



I. 

?p 

^3 

(/vi«Vj 

^ SCki o 

fm«vj 

1 

-O'HO 

i-M l.8o 2’53 

'-hOl -t*77 

-l‘7t 

t? 

-o-l^ 

O'fe? -4/*? 

//<• 

^?*7 

0*4 < Kc )*3 W' 
2*0 < K<.4*2?>*I* 

-0*37 

1‘^ 2.04 7 Hl 

-|.oo -|‘32 

-l 3 z_ 


-0*02 

I'l3 -3?'? 

KU. 

7‘?^7 

0*7 k.< 1-3 ^w::' 

2-0 ^m. < ^ ^4 

~0-3<f 

\'t Mo 2 

'*•{.00 —{.00 

— 1*00 

o 

-p'Oi 

1-20 -25 Z 

/27- 

74*0 

0*8^ te< l<4fvw' 

. «-l 
l.^ < te < T 

~D‘iQ 

2>^ 2'£>0 2 ‘ G>2 

-1.00 -0-77 

-Oil 

o 

-0 *oz 

-6i*30 

146 ■ 

67*7 

2l6. (e > |.0 

^0Z7 

2^2 2 8? 2 bl 

— j.oo -O'tl 

-0-hl 

o 

-0*02 

11-2 

164. 

64*? 

2t^. M. >M fv»:r' 

-O'Zi 

2-5^ 3'25 2-64 

-|.oo -0-44 

-P'44 

o 

-P-02 

I’W 30-4 

I?3. 

5f./ 

242. k > 1*3 fv«^‘ 

-0'2Z 

2-S 2 i,5 

-|.oo -0'33 

-0-33 

o 

-O’O/ 

1*0? 55T.3 

2o8. 

54^? 

262. h > l'3f^Jr' 

'^021 

3.0 3'to 2-fc7 

-1,00 -0*2? 

-D'Zg 

o 

-0-O« 

/'02 72*7 

225. 

S2-S 

27^. k>i'2 9v;:* 

-0-1$ 

3 2 42? 2 ?3 

-1*00 “0*23 

-0*25 

o 

-0-0/ 

0^1? /o8- 

2f>o. 

^7*3 

30?. k>»'iK* 

-on 

3 -5 4*44 2.?4 

— i'Oo — 0'l8 

-0-t% 

c> 

-0*01 

0-lfc ^34. 

2??. 

45'? 

334, k. > O.t^ f;' 

m 

3*? 5'Ot 24o 

— l*oO — 0*t5 

-0*/S" 

o 

-O.OJ 

0*^3 /72. 

325. 

43* S 

56?. k>0'7f>.T' 





TABLE A- 6 








Results from Method 3D 







M£TH0T> 

r 

V 

C-fwT') 

X 

y 

3P 

- 0'30 

2'0 

2 <■60 

2'6>Z 

3C 

O' 

2-2 

2'?0 

255 

3B 

O'Og 

4.0 

3'2 

0'«7 

3A 

~ 0'%3 

5'0 

4-0 

0'fc% 

^D 

“D‘7^ 

4.0 

4.0 

i«2^ 


±c. ^ 2 . ^3 

cmw; fWev) fMev) 


1*15 

-b'30 



216, 

l'»Z 

3'6? 

I5fc. 

67. 4 

223 


-27.0 

US. 

7?.l 

200. 

0'^6 

-18-2 

134. 

m.6 

212. 

0-‘^fc 


143, 

76'0 

224. 


TABLE A- 7 


Comparison of Results from Different Methods 


a 


PROJllCT B 

Bffects of Polarization on Neutron Star Structure 
B-1 Introduction 

Th? invisible components in the so called "single- line spectro- 
scopic binaries" in Hercules, Scorpius, etc., seen by the UHURU satel- 
lite are now commonly believed to be rotating neutron stars and in 
some cases may be even black hole revolving around the center of mass of 
the system. It is possible to estimate the mass of the neutron star 
component by an elaborate study of the intensity curves and of the 
spectral class of the visible components. Hence theoretical determinat- 
ions of masses, mpments of inertia and radii of stable neutron stars have 
become more important than ever. 

Macroscopic neutron star properties have been calculated during the 
past fifteen years tising equations of state derived from different pheno 
menological two-body interactions. One of the realistic effective inter- 
actions between two neutrons is a combination of a strong short-range re- 
pulsion and a long range attraction. It was first suggested by Migdal 
and later proved by Yang and Clark and others, that in a comparatively 

low-density degenerate neutron liquid, for which the interparticle spacing 

-13 

is large compared to the range of the repulsive forces (10 cm) , the 
attraction between pairs of neutrons of opposite spin and momentum would 
lead to the formation of a condensate and the appearance of superfluidity. 
However, our estimation of the pairing energy for "S" wave attraction 
is on the low side, since the enhancement of the attractive interaction 
between neutrons arising from the fact that they are embedded in a highly 
polarizable medium (i^he other neutrons) were not taken into account by us. 
(Yang and Clark), 
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B-2 Polarization Effect 

According to Pethick and Pines the additional terra coming from 
the polarizability of the medium is always attractive and is approximately 

- IVF^//fl ♦ F^) where V is the •'bare" interaction in the "S" wave 

' o‘ o 

channel and F^ is the Fermi liquid parameter which describes the spin- 
symmetric part of the interaction between two quasi-particles on the 
neutron Fermi surface. The net effective interaction, hence, has the 
simple form 

Since F® for the neutron liquid is negative and according to Pethick & 

o * 

Pines may be as negative as -0.7, such enhancement effects can be very 
important. An exact estimation of seems unattainable at present al- 
though such a calculation is desirable and necessary in order to under- 
stand exactly how the polarized medium ^.ffects the energy state and 

therefore the mass-energy density inside the neutron-star matter. 

- ^s 

B-3 An Approximate Calculation of F^ 

Applying the Landau technique of functional differentiation of the 
energy with respect to quasiparticle occupation numbers ^ ^ we estimate 
the Fermi liquid parameter F^ by summing over both the direct and exchange 
interactions between the interacting quasiparticle pair via the following 
integral 

pS Mo ^ \ VCo-> - VCk”>3 2.) 

^ -z ^ 

Where N = ra*k_/27^'fi^ is the familiar density of states at the Fermi surface, 
o f 

and iX is the volume of the system. 
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i’or the interaction between a quasiparticle pair, we choose the simple 
Yamaguchi potential^ 

(k'JVlk) = -(hVm) K g(k') g(k), (S-wave only) 
with gCk) = + k^)"^. For K = 0.18725 and 1,254 the singlet 

r.cattering length and effective range are a = -23.75 fm and r - 2.499 fitu 

S . . rrn. 

The Fermi liquid parameter is calculated for several densities. The 

results are listed in Table B-1 

'TUI r> 1 fof different densities 
Table B-1 o 


k^Cfra"^) 

0.(36 

0.60 

0.72 

0.96 

1.20 



-0.081 

-0.262 

-0,364 

-0.552 

-0.734 


0 





. .. . 



As can be seen from Table B-1, F^ depends on density quite strongly. As 
the density increases, the F^ becomes more and more negative. At the 
density k^ = 1.20 fm"^, is equal to -0.73, These results together with 
the prediction by Pethick and Finds give a very strong indication that the 
enhancement due to the polarization of the neutron medium may even be large 
enough to give rise to a major mass-energy density discontinuity due to 
neutron pairing effects . The existence of a concavity in the mass-energy 
density vs. the number density curve is sufficient to give a first-order 
phase- transition. Since there is no precise way of calculating the energy 
state for the neutron-etar matter including the polarization effect, it is cf 
interest to see how strong the enhancement will have to be in order to prodi.ice 
a first-order phase-transition. 
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B-4 Effect of Polarization on Condensation Energy and Superfluid-State Energy 

As has been mentioned in our original proposal we have all the niccb.saiv 

ingredients for the calculation of the normal state energy, and the 

condensatinn energy, £ , and also the superfluid state energy, £ (= £ — £ 

^ s n 

To simplify the calculation, we assume that the enhancement due to polaii- 
zation results in an increase in the attractive potential well-depth only. 

Then all we have to do in the calculations of the enhanced normal state energy, 
, and the enhanced condensation energy, £, , is to substitute A by 

AAq in the Ohmura potential 


V(12) = {SO 


r^ 2 < c 


S-wave only, ^ ^ 

enh enh , f enh 


I^A^ exp (-Kr), S-wave only, r^^^ ^ ( 8- ) 

Where c is the radius of the hard core. andP®”^ 

n ) 

-enh 

( =fcn — — ) are then calculated according to the procedures described in 

our proposal for the c = 0.4 fra Ohmura potential. The results of 5®”^, 

and ~ 1-0, 1.30, 1.45, and 1.50 are listed in Table B-2. 

(Note that a /? = 1.50 is corresponding to a = -0.33.] For the purpose of 
easier reference, we show in Figure B-1 the plots of the enhanced normal state 
energy per particle vs. for ( 3 - 1.0, 1.15, 1.30, 1.45, 1.50, 1.55, 1.60, 
and 1.80. In Figure B-2, the enhanced superfluid state energy per^ particle 
vs for 1.30, 1.45, and 1.50 are plotted. Some interesting sets 
of energies in Figure B-1 are plotted vs, the specific volume, iT, in Figure B-3. 
B-5 Discussion 

From Table B-2 and Figure B-2 we find that at ]3= 1.45, the superfluid state 
energy, turns negative at k^ = 0.5 fm"^, while at yS- 1.50, turns 
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negative at = 0.65 fin'^ resulting in a major mass-energy density dis- 

12 13 '3 

continuity around the density 10 -10 gm-cm . 

We conclude that the polarization effect indeed enhances the condensation 
energy (and the gap) and there is a tendency of the neutron-star matter to under- 
go a first-order phase- transition around the density of 10^^- 10^ gm-cm 
provided the effect is as strong as indicated (or stronger) . 
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Table B-2 Condensation energy, ^ , normal state energy per particle, and 
superfluid state energ^ per particle for 6 = 1 .0 Cunenhanced) , 
1.30, 1.45, and 1.50 (enhanced) 


k^(fm' 


n 


(MeV) 


(MeV) 


(MeV) 


1.0 

0.24 

0.48 

0.60 

0.72 

0.84 

1.08 

1.20 

0.606 

2.050 

2.960 

3.990 

5.150 

7.935 

9.655 

0.114 

0.294 

0.301 

0.220 

0.131 

0.012 

0.001 

0.492 

1.756 

2.659 

3.770 

5.019 

7.923 

9.654 

1.30 

0.36 

1.103 

0.436 

0.667 


0.48 

1.705 

0.773 

0.932 


0.60 

2.319 

1.120 

1.199 


0.72 

2.925 

1.030 

1.895 


0.84 

3.527 

0.765 

2.762 


0.96 

4.155 

0.504 

3.651 


1.08 

4.853 

0.352 

4.501 


1.20 

5.674 

0.107 

5.567 


1.32 

6.679 

0.021 

6.658 


1.45 

0.36 

0.48 

0,60 

0.72 

0.84 

0.96 

1.08 

1.20 

1.32 

1.028 

1.533 

1.996 

2.390 

2.719 

3.012 

3.315 

3.686 

4.189 

1.065 

1.601 

1.777 

1.685 

1.331 

0.933 

0.564 

0.277 

0.051 

-0.037 

-0.068 

0.219 

0.705 

1.388 

2.079 

2.751 

3.409 

4.138 

1.50 

0.36 

1.003 

1.270 

-0.267 


0.48 

1.476 

1.853 

-0.377 


0.60 

1.888 

2.070 

-0.182 


0.72 

2.212 

1.959 

0.253 


0.84 

2.450 

1.554 

0.896 


0.96 

2.632 

1.105 

1.527 


1.08 

2.803 

0.686 

2.117 


1.20 

3.023 

0.359 

2.664 


1.32 

3.359 

0.068 

3.-291 
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